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AB5TRACT 


An  investigation  of  using  the  response  of  a  circuit  at  selected  test 
frequencies  to  isolate  faulty  circuit  components  is  made.,  fc  procedure 
using  a  sensitivity  approach  for  intelligent  selection  of  test  frequencies 
is  developed.  The  developed  procedure  is  tested  and  the  results  are 
compared  with  results  using  conventional  procedures;..  The-  effect  of 
random,  within  tolerance  variations  of  nonfaulty  components:  on  the 
results  is  studied  for  both  conventional  and  developed  procedures . 
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I.  INTRODUCTION 

Maintenance  of  electronic  equipment  is  a  problem  that  increases  as 
equipment  becomes  more  complex.  With  ever  increasing  costs  (of  equip- 
ment) it  also  becomes  imperative  to  decrease  equipment  downtime  as 
much  as  possible.  To  this  end  a  fast,  reliable  method  of  locating 
circuit  failures  must  be  used. 

It  has  been  customary  to  isolate  failures  by  testing  individual  compo- 
nents. This  is  usually  done  by  isolating  the  fault  to  a  particular  stage, 
followed  by  subassembly  isolation  and  finally  component  isolation  by 
individual  testing. 

With  the  availability  and  flexibility  of  digital  computers,  a  pos- 
sible solution  to  the  problem  of  excessive  downtime  is  to  use  a  computer 
to  help  isolate  the  fault.  The  first  step  in  accomplishing  this  is  the 
development  of  a  testing  and  isolation  procedure  which  could  be  program- 
med for  the  digital  computer. 

A  recently  presented  procedure  has  been  examined  by  a  number  of 
individuals  [1],  [2],  [3],  and  [4].  The  procedure  was  described  by  Seshu 
and  Waxman  [1]  with  other  investigators  testing  and  making  applications  of 
the  procedure  with  few  modifications. 

The  basic  procedure  is  to  examine  the  response  of  a  circuit  under  test 
at  selected  frequencies.  Comparisons  of  the  results  with  a  list  of  "library 
results"  can  pinpoint  a  faulty  component. 

Most  investigators  have  used  test  frequencies  near  the  poles  and  zeros 
of  the  function  under  consideration.  Other  frequencies  have  been  chosen 
between  the  pole  and  zero  frequencies  [2]  and  [4].  In  this  study  an 
alternate  method  of  test  frequency  selection  was  developed  and  tested  and 
results  were  compared  with  results  using  other  methods. 
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The  procedures  for  signature  comparison  used  by  other  investigators 
have  not  been  explained  in  detail  in  most  reports.  However,  it  appears 
that  a  nearest  neighbor  approach  is  the  most  common  technique  used  to 
date.  Two  alternate  methods  were  examined  and  are  discussed  with 
conclusions  as  to  their  practicality. 

Previous  methods  of  fault  determination  after  signature  comparison 
have  also  not  been  documented.  Two  possible  methods  were  tried.  The 
first  method  used  computer  selection  for  each  network  function  and  manual 
correlation  for  fault  isolation.  The  second  method  was  completely  program- 
med for  digital  computation. 

As  far  as  can  be  determined  from  available  reports  no  investigator 
approached  the  problem  of  error  due  to  in-tolerance  variations  of  com- 
ponents. This  study  included  an  investigation  of  that  problem  with  both 
the  conventional  and  the  proposed  procedures. 
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II.     .PROCEDURAL  BACKGROUND 

A.  FUNDAMENTAL  CONCEPT 

All  components  of  a  given  network  play  a  part  in  the  formation  of  any 
of  the  network  functions.  The  sensitivity  of  any  given  network  function 
to  a  particular  component  is  generally  unique  to  that  component.  That  is, 
in  a  network  with  K  components  a  network  function  will  have  K  sensitivity 
functions,  one  for  each  component.  A  change  in  the  value  of  a  particular 
component  may  cause  a  decrease  in  the  magnitude  of  a  network  function  at 
one  frequency  and  an  increase  in  the  magnitude  at  a  different  frequency. 
There  may  in  fact  be  frequencies  at  which  the  function  sensitivity  is 
near  or  equal  to  zero  for  a  particular  component,  excluding  catastrophic 
failures. 

Specific  knowledge  about  this  behavior  can  be  used  to  detect  and 
isolate  component  value  changes.  If,  for  example,  it  is  known  that  the 
input  impedance  is  insensitive  at  low  frequencies  to  changes  in  a  specific 
component,  one  would  assume  that  component  was  not  faulty  if  the  input 
impedance  exhibits  a  significant  change  at  low  test  frequencies. 

B.  POLE-SHIFT  TECHNIQUE 

The  component  variation  effect  on  frequency  response  can  be  observed 
with  the  pole-shift  technique.  A  network  function  can  be  represented  by 
a  constant  and  a  unique  set  of  poles  and  zeros;  i.e., 


(s-Z^s-Zg)  •••  (s-Zm) 
F(S)  =  K  (s-Pl](s-p2)  ...  (s-pnj 


where 


F(s)  =  any  network  function 
K  =  constant 
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z.  =  i   zero  frequency 

p.  =  j   pole  frequency. 
J 

Every  component  plays  a  part  in  pole  and  zero  locations.  However,,  a 
few  components  may  dominate  the  location  of  a  given  pole  or  zero  frequency, 
Thus,  if  one  could  detect  the  shift  of  the  poles  and  zeros  of  a  gfven 
function,  it  may  be  possible  to  determine  component  variations.  A  simple 
example  can  best  illustrate  this  point. 
Example  1  -  Consider  the  network  shown  in  Figure  1. 
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Figure  1.  Example  1  -  Low  Pass  Filter 


v-i 

The  network  function  j—  (s)  is  given  by 


s'  +  lis  +  11 
s2  +  10s  +  1 


-  (s  +  9.887)(s  +  1.113) 
"  (s  +  9.899)(s  +  0.101)  ' 

The  pole  and  zero  shifts  were  determined  by  varying  each  component 
value  by  ten  percent  and  observing  the  corresponding  critical  frequency 
shifts.  The  results: 
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Location 

(radians 

/sec) 

Percent  Shift 

Circuit 
Condition 

pl 

^2 

zl 

z2 

h 

P2 

ZT 

z2 

Nominal 

.101 

9.899 

1.113 

9.887 

0 

0 

Q 

0 

R,   variation 

.101 

9.899 

1.020 

9.889 

0 

0 

JG2 

.02 

C,   variation 

.092 

9.908 

1.010 

9.899 

8.9 

.09 

9..Z5 

.12 

L,   variation 

.101 

8.990 

1.114 

8.977 

0 

9.2 

.as 

9.2 

R2  variation 

.092 

10.908 

1.101 

10.899 

8.9 

10.2 

T.07 

10.2 

TABLE  II-l.  Percent  Pole  and  Zero  Shifts 

It  is  immediately  obvious  that  a  change  in  R?  or  C,  will  cause  a  greater 
shift  in  z.  than  will  a  corresponding  change  in  L,  or  R, .  Also,  p?  and  z? 
are  most  affected  by  L-.  and  Rp  and  p,  is  affected  only  by  C,  and  R„. 

C.  CHANGES  IN  RESPONSE  FROM  POLE  SHIFT 

In  practice  the  exact  location  of  the  critical  frequencies  can  usually 
not  be  determined  directly  from  the  circuit.  However,  this  does  not  rule 
out  the  pole  shift  technique.  Any  shift  in  critical  frequency  causes  a 
change  in  the  magnitude  of  the  given  function  at  all  frequencies  greater 
than  the  shifted  one.  This  can  be  seen  from  the  simple  Bode  representation 
of  a  pole: 


|G(f) 


Figure  2.  Bode  Plot  of  Simple  Pole 
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If  f  is  changed  to  f'  (indicated  by  dashed  line),  the  function 
magnitude,  at  any  frequency,  f  ,  greater  than  f  ,  is  changed  (increased 
in  this  case).  The  magnitude  is  not  affected  at  frequencies  less  than 


V 


Now,  from  the  network  in  Figure  1  we  obtain: 
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10 


*1   s2  +  lis  +  11    (s  +  9.887)(s  +  1.113) 

Figure  3a  is  the  Bode  plot  for  the  function  Vp/V-, .  Figure  3b  and 
3c  show  the  effect  of  increasing  R,  and  of  increasing  L, ,  respectively, 
The  nominal  response  is  shown  as  a  dashed  line. 


H(w) 


10 

0 

-10 


db 


.01 


p1  =  .099 


p2  =  1.57  Hz 


.1  1. 

Figure  3a.     V2/V]   Bode  Plot, 


10. 


.1  1.       f 

Figure  3b.  V"2/V, ,  R,  Increased. 


.1  1.      f 

Figure  3c.  V^/V-j ,  L1  Increased 
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It  can  be  seen  that  an  increase  in  R-,  causes  a  decrease  in  the 
magnitude  of  the  response  between  p,  and  p2-  A  lesser  decrease  is  seen 
at  frequencies  greater  than  p2-  The  response  between  p,  and  p?  is 
increased  for  an  increase  in  L,  while  the  response  above  p«  is  decreased 
Thus,  if  one  examines  the  response  at  a  frequency  between  p,  and  p2  and 
at  a  frequency  greater  than  p2  it  will  be  a  reasonably  simple  task  to 
distinguish  between  an  increase  in  R,  or  L, .  This  example  is  an  over- 
simplification of  the  procedure,  but  the  application  of  which  can  be 
appreciated. 
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III.  GENERAL  PROCEDURE 

A.  BASIC  PROCEDURE 

The  basic  procedures  outlined  by  Seshu  and  Waxman  [1]  was  followed 
with  modifications  as  explained  herein.  Two  network  functions  were  used 
concurrently.  Work  done  by  Maenpaa,  Stehman,  and  Stahl  [5]  indicated 
that  the  redundancy  from  using  two  functions  was  more  than  adequately 
compensated  for  by  the  increased  fault  detection  capability. 

The  symbolic  network  functions  were  calculated  by  hand.  It  should 
be  pointed  out  that  computer  programs  exist  which  calculate  symbolic  net- 
work functions.  These  programs  should  be  used  for  more  complex  circuits. 

The  functions  were  numerically  evaluated  and  factored  to  determine 
the  poles  and  zeros.  Test  frequencies  were  then  chosen.  Previous 
procedures  generally  selected  test  frequencies  based  on  the  pole  and 
zero  locations.  As  a  test  of  the  procedure  several  test  frequencies  were 
chosen  between  the  critical  frequencies  as  well  as  0  hertz  and  one 
frequency  above  the  highest  critical  frequency.  The  actual  selection 
was  programmed  for  computer  selection  and  is  described  in  Section  VII. 

It  was  felt  that  the  procedure  of  selecting  test  frequencies  based 
on  poles  and  zeros  was,  although  intuitively  satisfying,  a  somewhat 
illogical  procedure.  The  concern  was  not  with  the  function  itself  but 
with  the  behavior  of  that  function  with  a  parametric  change.  Thus, 
sensitivity  functions  were  calculated  for  each  element.  Based  on  these 
functions,  six  test  frequencies  were  chosen.  This  procedure  is  discussed 
in  Section  IV. 

After  selection  of  test  frequencies,  a  worst  case  response,  for 
component  values  within  10  percent  tolerance  limits,  was  computed  for 
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each  test  frequency.  This  defined  a  nominal  range  of  performance.  This 
is  not  necessarily  the  best  method  for  definition  of  nominal  range.  This 
range  might  be  specified  for  a  particular  circuit.  However,  in  the 
absence  of  such  specifications  the  above  procedure  was  followed. 

Each  component  was  then  allowed  to  take  on  several  discrete  values 
out  of  tolerance  limits.  All  other  components  were  held  fixed  at  their 
nominal  value.  The  response  at  each  test  frequency  was  then  calculated 
and  catalogued  for  the  entire  collection  of  sets  of  component  values. 

The  catalogued  response  values  were  then  quantized  into  nine  levels 
and  labeled  one  through  nine.  Quantization  level  five  was  reserved  for 
response  within  nominal  limits.  The  setting  of  the  levels  is  explained 
in  Part  B.  The  quantized  responses  for  a  given  set  of  component  values 
is  defined  as  the  signature  for  the  given  circuit  configuration. 

The  signatures  were  computed  for  the  following  cases: 

value  =  zero 
+  20% 
+  50% 
100  times  nominal . 

The  entire  set  of  signatures  for  a  given  function  is  defined  as  the 

library  for  that  function.  There  is  no  theoretical  limit  on  the  size 

of  the  library.  However,  the  inclusion  of  many  more  circuit  configurations 

would  probably  only  complicate  the  matching  procedure  and  not  increase 

detection  capability. 

Using  the  primary  test  circuit  the  input  impedance  and  voltage  gain 
at  each  test  frequency  was  measured.  Then,  after  quantization,  the  signa- 
tures measured  were  compared  with  the  library  signatures  and  a  "closest 
match"  was  made.  Thus,  the  faulty  component  was  identified. 

In  order  to  test  the  procedure,  several  error  conditions  were  set  and 
the  response  calculated.  In  order  to  more  closely  approximate  actual 
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circuit  conditions,  all  component  values  were  allowed  to  vary  randomly 
(uniform  distribution)  within  10  percent  tolerance  limits.  Calculations 
were  made  with  and  without  the  random  variations  to  determine  how  much 
effect  this  would  have  on  the  results. 

The  results  of  all  tests  are  presented  in  Section  VI. 

B.  SETTING  OF  QUANTIZATION  LEVELS 

.  With  the  exception  of  selection  of  test  frequencies  the  most  critical 
step  in  the  procedure  was  the  setting  of  quantization  levels.  Two  methods 
were  used  successfully.  One  was  programmed  as  part  of  the  main  analysis 
program.  The  second  method,  although  not  programmed  is  possibly  adaptable 
to  computer  selection.  Both  methods  have  nine  levels  of  response  with 
level  five  reserved  for  the  nominal  range. 

The  nominal  range  was  selected  by  taking  the  range  between  the 
.maximum  positive  and  maximum  negative  deviation  from  normal  for  any 
clement  value  at  +10%  from  normal.  The  programmed  method  then  set  the 
other  eight  levels  as  shown  in  Table  III-l  in  which  F  is  a  multiplying 
factor  used  to  force  the  possible  response  range  to  cover  all  quantum 
levels  with  levels  one  and  nine  assigned  only  for  extreme  cases. 


Quantum  Level 

Normalized  Response 

1 

(oo,  1+13  FN] 

2 

(1   +  13  FN,  1   +  7  FN] 

3 

(1   +  7  FN,  1   +  4  FN] 

4 

(1   +  4  FN,  1   +  N] 

5 

(1   +  N,  1    -  M) 

6 

[1   -  M,  1    -  4'FM) 

7 

[1   -  4  FM,  1    -  7  FM) 

8 

[1   -  7  FM,  1    -  13  FM) 

9 

[1   -  13  FM,  0) 

TABLE  III-l.  Quantization  levels. 
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The  second  method  used  a  more  inituitive  approach.  The  maximum 
deviation  levels,  1  and  9,  were  selected  based  on  the  maximum  deviation 
from  normal  under  any  condition.  The  mid  levels  were  then  selected  so 
that  the  catalogued  or  library  signatures  would  cover  all  quantum  levels 
as  uniformly  as  possible. 

The  manual  approach  is  acceptable  for  small  circuits  but  gets 
unmanageable  quickly  as  circuit  size  increases.  Parts  of  the  procedure 
could  be  easily  programmed.  However,  a  suitable  algorithm  would  have  to 
be  developed  prior  to  complete  programming. 

C.  TWO  COMPONENT  VARIATIONS 

It  cannot  be  assumed  that  only  one  component  at  a  time  will  undergo 
changes  in  value.  In  order  to  investigate  the  possibility  of  two  compo- 
nent values  changing  simultaneously,  two  possibilities  were  considered. 
If  one  component  masks  function  changes  due  to  variations  in  another 
component  (i.e.,  no  change  in  function  magnitude),  the  entire  procedure 
will  not  work.  However,  if  one  is  primarily  concerned  with  circuit 
performance,  the  failure  to  detect  component  variations  is  of  no  concern 
if  there  is  no  degradation  in  circuit  performance. 

The  second  situation  considered  was  one  in  which  there  was  an  obvious 
change  in  circuit  response  due  to  two  components  varying.  The  original 
procedure  made  allowances  for  varying  only  one  component  at  a  time.  How- 
ever, a  slight  modification  allowed  for  the  generation  of  a  set  of  library 
signatures  for  two  component  variations.  The  set  is  presented  in  the 
computer  output  section. 

It  soon  became  apparent  that  in  order  to  use  the  two  component 
variation  library,  one  needed  apriori  knowledge  of  the  number  of  faulty 
components.  Since  this  is  not  usually  available,  the  usefulness  of  two 
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component  signatures  is  doubtful.  Additionally,  an  examination  of  the 
library  signatures  revealed  that  one  of  the  two  components  varied  usually 
dominates  the  response.  Thus,  if  two  components  have  failed  the  single 
component  library  would  probably  isolate  the  dominant  component.  After 
replacement  of  that  component  further  testing  could  reveal  the  second 
faulty  component.  Consequently  no  further  investigation  of  multi- 
component  variation  was  undertaken. 
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IV.  PRIMARY  TEST  CIRCUIT 

The  primary  circuit  used  throughout  the  investigation  was  a  six- 
element  low-pass  filter.  (See  Figure  4.)  The  unsealed  frequency 
response,  represented  by  the  voltage  ratio  transfer  function,.  V2/Vy,, 
is  shown  in  Figure  5  and  the  unsealed  input  impedance  as  a  functiorr 
of  frequency  is  shown  in  Figure  6.  On  Figures  5  and  6  the  poles  are 
represented  by  X  and  the  zeros  by  0. 
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Figure  4.  Primary  Test  Circuit. 


The  component  values  were  magnitude  scaled  by  100  and  frequency 

5 
scaled  by  10  .  This  yielded  component  values  of  approximately  the: 

same  magnitude.  Ideally,  the  scaled  cutoff  frequency  shouTd  be  T  Hz. 

However,  with  this  scaling  f  =  .02  Hz  and  no  difficulty  was  encounter- 


ed 


The  scaled  network  functions  of  interest  are: 

.1634 


'2  = 
^l"      s4+17.67s3+6.462s2+2.601s+.2124 


.1634 

(s+17.31)(s+.0997)(s+1.325+j3.25)(s+K325-j3.25) 
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Figure  5.  Vp/V-, ,  Primary  Test  Circuit 


.01  .1  1. 

Figure  6.  V,/!-,,  Primary  Test  Circuit 
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V1   _     s4+17.67s3+6.462s2+2.601s-K2124 
h         5.882s3+.147s2+.817s+.01634 


s+17.31)(s+.0997)(s+1.325+j3.25)(s-H.325-j3.25) 
s+.02)(s+.0025+j.373)(s+.0025-j.373) 


The  poles  and  zeros  in  Hz  are  listed  in  Table  IV-1  for  convenience, 

V1!   .  V1! 

Zeros*  (Hz)  Poles  (Hz) 

.0518  (Second  Order)    .0594  (Second  Order) 

.0159  .00319 

2.76 

*  Zeros  of  V^/I,  are  the  poles  of  VyV-j . 

TABLE  IV-1.  Pole  and  Zero  Magnitudes. 
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V.  SELECTION  OF  TEST  FREQUENCIES 

A.  CONVENTIONAL  METHOD 

The  conventional  procedure  was  followed  with  regard  to  selection  of 
test  frequencies  in  order  to  test  the  method.  It  soon  became  apparent 
that  the  number  of  test  frequencies  was  too  large  to  handle  with  ease. 
Thus,  only  frequencies  between  the  pole  and  zero  frequencies  were 
selected  as  well  as  0  Hertz  and  one  above  the  maximum  criticaT  frequency. 
The  computer  program  written  for  implementation  of  the  procedure  (see 
Program  1)  did  the  frequency  selection  automatically.  The  test  fre- 
quencies were  selected  by  taking   .,,    . 

.  f1  '  -  f1 
f  =  f 1  +  c  «  c 
Ttest   Tc  -    2 

for  ewery  f    where  f  =  i   critical  frequency.  If 


,    fi+1  -  f1' 


:i+l 


fi+2       fi+l 
c  c 


one  of  the  frequencies  was  eliminated.  Test  frequencies  were  computed 
for  the  input  impedance  because  the  critical  frequencies  for  V^/V,  are 
included  in  the  critical  frequencies  for  V,/I,.  One  test  frequency  was 
eliminated  due  to  the  above  inequality.  This  resulted  in  the  selection 
of  the  following  seven  test  frequencies  for  a  fourth  order  filter: 


0.0 

Radians 

.001 

Radians 

.00498 

Radians 

.196 

Radians 

.5488 

Radians 

8.49 

Radians 

25.8 

Radians 

TABLE  V-l .  Conventional  Method  Test  Frequencies 
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B.  SENSITIVITY  APPROACH  TO  FREQUENCY  SELECTION 

Inspection  of  the  signatures  generated  using  the  conventional  method 
showed  that  redundancy  existed  to  the  extent  that  the  method  was  some- 
what inefficient.  Two  test  frequencies,  D.C.  and  the  lowest  nonzero  test 
frequency,  gave  nearly  identical  results.  This  indicated  one  of  the  two 
frequencies  v/as  probably  not  required.  It  was  also  felt  that  unless 
further  investigation  was  undertaken  it  would  be  possible  to  ignore  a 
range  of  frequencies  which  could  be  extremely  helpful. 

This  investigation  took  the  form  of  calculating  and  plotting  functions 
which  defined  the  sensitivity  of  V«/V1  and  I-i/V,  to  variations  in  each 
component.  These  plots  are  presented  as  Figures  7  through  18.  Figure  19 
is  a  superposition  of  the  plots  of  the  sensitivity  function  for  1-,/V,  and 
Figure  20  is  the  same  for  Vp/V, .  Ideally  V,/I  should  have  been  used. 
However,  the  computer  program  used  in  this  portion  of  the  study  did  not 
allow  for  using  V-,/1-,.  It  was  felt  that  no  loss  of  information  would 
result. 

The  superposition  plots  Figures  19  and  20  are  very  revealing.  The 
first  observation  is  that  Vp/V,  is  not  sensitive  to  changes  in  R-,  and  R2 
at  high  frequencies  nor  to  C's  or  L's  at  low  frequencies.  While  this 
information  is  not  new  the  sensitivity  functions  assist  in  determining 
when  C's  and  L's  have  little  effect.  This  can  be  used  to  eliminate  the 
redundancy  observed  using  the  conventional  method  of  frequency  selection. 

The  mid-frequency  range  is  the  most  interesting  range  for  the  selection 
of  test  frequencies.  For  example,  near  .02  Hz  V?/V,  is  very  sensitive  to 
changes  in  R,  and  C2»  less  sensitive  to  changes  in  L2»  even  less  sensitive 


*  The  program  NASAP-70  v/as  used  to  calculate  the  various  sensitivity 
functions. 
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to  changes  in  R2  and  C,  and  virtually  insensitive  to  changes  in  L, . 
There  is  a  peaking  of  the  sensitivity  functions  for  C, ,  C?,  and  L?  near 
.06  Hz.  At  this  same  frequency  the  sensitivity  functions  for  R-, ,  R?, 
and  L,  are  near  zero. 

The  I/,  and  V-,  plots  show  that  at  near  .3  Hz  only  the  sensitivity 
function  for  R-,  is  high.  The  functions  for  L^  and  C-,  are  low  and  the 
other  functions  are  zero.  At  3  Hz  T-i/V-,  is  about  equally  sensitive  to 
changes  in  L,  and  R-.  while  the  other  functions  are  very  close  to  zero. 
It  can  also  be  seen  that  the  frequencies  discussed  in  connection  with 
Vo/V-j  are  reasonably  good  choices  for  I-i/V-i. 

Based  on  the  above  observations  the  following  six  frequencies  were 
selected  as  test  frequencies: 

Number  Frequency 

1  .00048  Hz 

2  .027  Hz 

3  .06  Hz 

4  .299  Hz 

5  2.99  Hz 

6  49.97  Hz 

TABLE  V-2.  Sensitivity  Function  Test  Frequencies. 

After  some  preliminary  work  it  was  observed  that  the  response  at  .06 
Hz  was  not  well  behaved.  Since  four  sensitivity  functions  peaked  for 
1-,/V-j,  random  variation  within  tolerance  limits  of  these  components  were 
being  observed  as  significant  changes  in  response.  In  order  to  overcome 
this  problem,  test  frequency  three  was  changed  to  .055  Hz.  The  new  fre- 
quency was  off  the  peak  enough  to  be  useable  and  still  yield  the  desired 
information. 


26 


4. 


Fs(w) 


3. 


2. 


1. 


0. 


.001 


.01 


.1 


10. 


w 


Figure  7.  Sensitivity  of  Vp/V,  to 
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Figure  8.  Sensitivity  of  V2/V1  to  LI 
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Figure  9.  Sensitivity  of  V^/V,  to  CI 
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Figure  10.  Sensitivity  of  Vg/V,  to  R2 
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Figure  11.     Sensitivity  of  vyv,   to  L2 
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Figure  12.     Sensitivity  of  vyV-j   to  C2. 
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Figure  13.     Sensitivity  of  I,/V,   to  Rl . 
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Figure  14.     Sensitivity  of  1,/V-,   to  LI. 
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Figure  15.     Sensitivity  of  I,/V,   to  CI 
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Figure  16.     Sensitivity  of  I,/V|   to  R2- 
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Figure  17.     Sensitivity  of  I,/V,   to  L2. 
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Figure  18.     Sensitivity  of  I,/Vj   to  C2, 
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Figure  19.  I, /V-.  Composite  Sensitivity  Functions 
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Figure  20.  Vp/V-,  Composite  Sensitivity  Functions 
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C.  COMPARISON  OF  METHODS 

The  comparison  of  the  results  using  the  two  methods  of  frequency 
selection  is  included  in  Section  VI.  Some  comments  as  to  applicability 
should  be  made  at  this  point. 

If  one  computes  sensitivity  functions  by  using  the  derivative:  of 
the  network  function  with  respect  to  the  element  under  consideration  it 
becomes 


fs(s)  = 

d         p(s) 
dek      q(s) 

q(s)     Mi).p(s)Mil 
k                        k 

[q(s)]2 

where: 

fs(s)  =  sensitivity  function 

d/de.  =  derivative  with  respect  to  element  k 

Pr   =  ratio  of  polynomials  representing  the  network  function 

Thus,  the  poles  of  the  sensitivity  function  are  the  same  as  the.  poles  of 
the  original  function.  Consequently,  if  one  selects  test  frequencies  on 
the  basis  of  function  poles  the  selection  is  linked  to  the  sensitivity 
functions.  Generally,  however,  the  zeros  of  the  sensitivity  function  are 
unique  and  are  thus  ignored  by  any  selection  procedure  based  on  the  poles 
and  zeros  of  the  original  function. 

From  the  preceding  paragraph,  some  of  the  advantages  of  using  sensi- 
tivity functions  for  test  frequency  selection  become  apparent..  However, 
that  method  has  the  disadvantage  of  requiring  the  calculation  orf  sensi- 
tivity functions  for  each  element. 
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With  the  small  circuit  used  as  an  example  the  value  of  sensitivity 
functions  are  demonstrated.  With  more  complex  circuits  it  appears  more 
probable  that  important  test  frequencies  might  be  overlooked  if  the  sen- 
sitivity functions  are  not  used. 


42 


VI.  RESULTS  WITH  PRIMARY  CIRCUIT 

A.  RESULTS  USING  BASIC  PROCEDURE 

Using  the  procedures  outlined  in  Section  III-A  and  V-A  a  set  of 
signatures  was  generated  for  the  six  selected  circuit  conditions. 
These  signatures  are  presented  in  Figure  25.  A  second  set  of  signatures 
was  then  generated  for  the  same  circuit  conditions  with  random  variation 
allowed  for  the  circuit  elements.  The  second  set  of  signatures  was  input 
to  the  selection  program  to  determine  what  errors  resulted  from  the  random 
variations. 

Table  VI-1  lists  the  input  condition  and  the  program  output.  The 
numerical  result  was  determined  as  follows: 

1  =  Selection  was  correct  condition 

2  =  Selection  was  correct  element 

3  =  One  of  multiple  selections  was  the  correct  condition 

4  =  One  of  multiple  selections  was  correct  element 

5  =  Incorrect  selection. 

By  correlating  the  results  from  Vp/V,  and  V./L  a  "hit",  "miss",  or 
"possible  hit"  rating  was  assigned  for  each  test  signature  as  follows. 
A  hit  was  assigned  if  both  V2/V,  and  V,/I,  yielded  a  1  or  2  or  if  one 
was  1  or  2  and  the  other  was  3  or  4.  A  3-5,  4-5,  or  5-5  combination  was 
called  a  miss.  All  other  combinations  were, called  possible  hits.  For 
some  cases  deviations  were  warranted.  For  example,  a  3-4  selection  could 
have  been  called  a  miss  if  a  component  other  than  the  given  component 
was  clearly  indicated  as  being  faulty. 

The  result  was  21  hits,  9  possible  hits  and  6  misses  for  58%  hits  and 
83%  hits  or  possible  hits.  Since  only  the  circuit  conditions  used  for 
signature  generation  were  used  to  generate  the  test  signatures  100% 


43 


Circuit 

Program 

Numerical 

Program 

Numerical 

Condition 

Con 

dition 

Result 

Condition 

Result 

Rl :  Short 

Rl: 

Short 

1 

Rl: 

Short 

1 

-50% 

Rl: 

-50% 

1 

L2: 

Open 

5 

-20% 

Rl: 

Open 

2 

L2: 

+50% 

5 

+20% 

L2: 

+20% 

5 

Rl: 

+20% 

1 

+50% 

Rl: 

+20% 

2 

Rl: 

+20% 

2 

Open 

Rl: 

Open 

1 

L2: 
Rl: 

Short 
Open 

3 

LI:  Short 

LI: 

Short 

1 

LI 

-50% 

2 

-50% 

LI: 

-50% 

1 

LI 

-50% 

1 

-20% 

LI: 
R2: 

-20% 
+20% 

3 

LI 

-20% 

1 

+20% 

LI: 

+20% 

1 

LI 

-20% 

2 

+50% 

LI: 

+20% 

2 

LI 

'  -20% 

2 

Open 

LI: 

Open 

1 

R2 

+20% 

5 

CI :  Open 

CI: 

Open 

1 

LI. 

+50% 

5 

-50% 

CI: 
C2: 

+20% 
-50% 

4 

CI: 

-50% 

1 

-20% 

C2: 

+20% 

5 

Rl 

-20% 

5 

+20% 

LI: 
C2: 

+20% 
+20% 

5 

LI 

-50% 

5 

+50% 

C2: 

+20% 

5 

CI. 
C2: 

+20% 
+20% 

3 

Short 

CI: 

Short 

1 

CI 

Short 

1 

R2:  Short 

R2: 

-50% 

2 

R2 

Short 

1 

-50% 

R2: 

-50% 

1 

R2 
R2 

-50% 
•  -20% 

2 

-20% 

R2: 
L2: 
C2: 

-20% 
+20% 
+20% 

3 

Rl 

:  -20% 

5 

+20% 

C2: 

+20% 

5 

CI 
C2 

+20% 
+20% 

5 

+50% 

R2: 

+20% 

2 

R2 

:  +20% 

2 

Open 

R2: 

Open 

1 

CI 

:  Open 

5 

TABLE  VI 

-1 .  Basic 

Procedure 

>  Results. 
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L2:  Short 

L2:  Short 

-50% 

L2:  -50% 

-20% 

L2:  -50% 

LI:  +20% 

+20% 

L2:  +20% 

C2:  +20% 

+50% 

L2:  +20% 

Open 

L2:  -20% 

C2:  Open 

CI:  +20% 

C2:  -50% 

-50% 

CI:  -20% 

-20% 

Nominal 

+20% 

C2:  +20% 

+50% 

Nominal 

Short 

C2:  Short 

R2:  Short 

L2:  Short 

Rl :  Open 

L2:  -50% 

L2:  -20% 

L2:  +20% 


L2:  Open 
L2:  Open 


R2:  +20% 

CI:  +50% 

CI:  -20% 

C2:  -20% 

L2:  +20% 

C2:  +20% 

C2:  +20% 

CI:  +20% 

C2:  Short 


4 
5 
3 
4 


TABLE  VI-1  (Continued) 
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hits  or  possible  hits  could  have  been  accomplished.  Thus,  using  the 
conventional  procedure,  within  tolerance  variations  introduced  17%  to 
42%  errors. 

B.  RESULTS  USING  SENSITIVITY  FUNCTIONS 

Using  the  frequencies  determined  from  the  sensitivity  functions  a 
set  of  signatures  was  generated  for  the  six  circuit  conditions.  In  order 
to  check  the  effect  of  variations  within  tolerance  for  this  method  the 
same  procedure  as  used  in  VI-A  was  followed.  The  result  was  50  hits,  16 
possible  hits  and  6  misses.  The  6  misses  included  3  due  to  the  test 
signature  showing  normal  response  from  the  circuit.  These  results  showed 
an  error  rate  of  11%  -  31%. 

Eight  additional  circuit  conditions  were  simulated  for  each  element 
to  see  what  could  be  expected  from  a  less  artificial  situation.  Table 
VI-2. lists  all  the  conditions  set  and  the  numerical  result.  The  numerical 
result  was  assigned  as  follows: 

1  -  closest  library  signature  selected 

2  -  right  element  selected 

3  -  one  of  multiple  signature  closest  to  given  condition 

4  -  one  of  multiple  signatures  right  element 

5  -  incorrect  selection. 

Hits,  misses  and  possible  hits  were  also  assigned  as  before.  The 
result  was  115  hits,  43  possible  hits  and  10  misses.  Or,  68%  hits  and 
94%  hits  or  possible  hits. 

C.  USING  SUBROUTINE  DET  FOR  CORRELATION 

The  same  set  of  test  signatures  used  in  VI-B  was  used  to  test  the 

programmed  correlation  and  selection.  The  input  and  results  are  listed 
in  Table  VI -3.  The  results  are  coded  as  follows: 
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Circuit 
Condition 

Rl :  Short 
-90% 
-70% 


Run  1 


Run  2 


-40% 


+20% 

+30% 

+50% 

+60% 

+90% 

x3 

X75 

Open 


LI :  Short 


-50% 

-40% 

-20% 

+20% 

+30% 

+50% 

+60% 

+90% 

X3 

X75 

Open 


Vvi 


V1! 

1 
1 
1 
1 
1 
5 
3 
1 
2 
1 
•1 
1 
1 
1 


Vvi 


V1! 


R2:  Short 
-90% 
-70% 
-50% 
-40% 
-20% 
+20% 
+30% 
+50% 
+60% 
+90% 
Open 


3 
3 
3 


TABLE  VI-2.  Results  Using  Sensitivity  Functions 
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CI :  Open 

1 

1 

1 

-90% 

5 

1 

5 

-70% 

5 

1 

5 

-50% 

3 

1 

5 

-40% 

3 

1 

3 

-20% 

1 

3 

5 

•+20% 

4 

5 

5 

+30% 

3 

1 

5 

+50% 

5 

1 

5 

+60% 

3 

1 

5 

+90% 

4 

1 

3 

X3 

5 

5 

5 

4 

X75 

1 

1 

1 

Short 

1 

1 

1 

L2:  Short 

2 

1 

2 

1 

-90% 

1 

5 

1 

5 

-70% 

1 

1 

1 

5 

-50% 

1 

4 

2 

3 

-40% 

1 

3 

1 

3 

-20% 

2 

5 

3 

3 

+20% 

3 

3 

4 

3 

+30% 

5 

3 

5 

3 

+50% 

5 

1 

3 

1 

+60% 

5 

1 

5 

1 

+90% 

5 

1 

5 

1 

X3 

5 

1 

5 

1 

X75 

5 

1 

5 

1 

Open 

1 

1 

5 

1 

C2:  Open 

1 

1 

1 

1 

-90% 

1 

1 

5 

1 

•  -70% 

1 

1 

5 

1 

-50% 

3 

1 

4 

1 

-40% 

3 

1 

3 

3 

-20% 

5 

3 

5 

3 

+20% 

5 

3 

5 

3 

+30% 

3 

1 

5 

1 

+50% 

5 

1 

5 

1 

+60% 

3 

1 

3 

1 

+90% 

3 

1 

5 

1 

X3 

5 

1 

5 

1 

X75 

3 

3 

3 

3 

Short 

3 

3 

5 

3 

TABLE  VI-2   (Continued) 
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1  -  correct  primary  selection 

2  -  one  of  two  primary  selections  count 

3  -  one  of  two  secondary  (no  primary  selections) 

4  -  incorrect  primary  selection. 

The  results  were  61  rated  1,  12  rated  2,  21  rated  3  and  2  rated  4. 
Or,  63%  completely  correct  and  98%  with  one  of  two  selections  correct. 
This  is  roughly  the  same  result  using  manual  correlation.  The  algorithm 
used  for  determining  primary  and  secondary  selection  is  outlined  in 
Section  VII-C. 
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Circuit  Condition 


Vvi 


V1! 


Rl 

-90% 


+30% 

+60% 

+90% 

X3 

X75 


LI 

-90% 


-60% 

+30% 

+60% 

+90% 

X3 

X75 


CI 

-90% 


-60% 

+30% 

+60% 

+90% 

X3 

X75 


3 
1 
3 
3 
3 
1 
4 
1 


3 
1 
3 
2 
2 
1 
4 
1 


R2 

-90% 
-70% 


+30% 

+60% 

+90% 

X3 

X75 


TABLE  VI-3.  Computer  Correlation  Results 
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L2 

-90°/  3  3 


3  1 

-60%  1  1 

+30%  2  2 

+60%  2  2 

+90%  3  3 

X3  3  3 

X75  3  3 


C2 

-90%  3  3 

-70%  2  2 

-60%  3  3 

+30%  1  1 

+60%  1  1 

+90%  1  1 

X3  3  3 

X75  2  2 

TABLE  VI-3  (continued) 
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VII.  COMPUTER  PROGRAMS 

Several  computer  programs  were  written  in  Fortran  IV  to  accomplish 
much  of  the  mathematical  work  associated  with  this  study.  Four  of  the 
programs  are  included  herein  with  detailed  descriptions.  Each  of  these 
programs  were  used  with  several  modifications.  However,  only  the  basic 
program  is  presented  in  each  instance.  All  listings  are  in  the  designated 
section. 

A.  PRIMARY  ANALYSIS  PROGRAM 

The  primary  analysis  program  contains  the  main  program  and  three 
subroutines.  A  fourth  subroutine,  PRQD,  is  used  for  root  finding  and 
is  not  included  in  the  listing. 

The  program  calculates  the  magnitude  of  both  network  functions  as  a 
function  of  the  complex  frequency,  s,  finds  the  roots  of  the  polynomials 
involved  and  selects  test  frequencies. 

Each  element  in  turn  is  then  set  to  each  out-of-tolerance  value  as 
selected  previously,  the  magnitude  of  each  function  at  each  test  fre- 
quency is  calculated  and  all  responses  are  quantized  into  signatures. 

The  output  consists  of  both  network  functions,  the  test  frequencies, 
the  nominal  response  at  each  test  frequency  and  the  list  of  library 
signatures. 

The  listing  is  the  program  which  calculates  the  test  frequencies  as 
explained  in  Section  IV-A.  To  modify  the  program  for  the  selection  from 
sensitivity  functions  the  portion  of  the  program  between  indicated  com- 
ment cards  was  replaced  with  a  read  statement  and  the  test  frequencies 
were  used  as  the  input. 

Figure  21  is  the  flow  graph  for  the  main  analysis  program. 

52 


NO 


READ 
QUANTI- 
ZATION 
FACTOR 


SET  NEW 
COMPONENT 
VALUES 


CALCULATE 
AND 


CALCULATE 

FUNCTION 

MAGNITUDE 


OUTPUT 
NOMINAL 
RESPONSE  AND 
TEST  FREQS. 


CALCULATE 

FUNCTION 

MAGNITUDES 


YES 


FIND  ROOTS 
AND  SELECT 
REST  FREQS 


YES 


CALCULATE 
NOMINAL 
RESPONSE 
RANGE 


QUANTIZE 
AND  FORM 
SIGNATURES 


OUTPUT 
SIGNA- 
TURES 


Figure  21.  Flowgraph  of  Primary  Analysis  Program 
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B.  TWO  COMPONENT  VARIATION  PROGRAM 

The  main  analysis  program  was  modified  to  generate  signatures  for 
simultaneous  two-component  failure.  Two  different  modifications  were 
made,  one  generated  signatures  for  only  short  and  open  element  failures. 
The  other  allowed  for  short,  -50%,  +50%  and  open  failures. 

Since  the  programs  were  only  modifications  to  the  main  analysis 
program  flow  data  are  not  included  and  the  listings  are  omitted. 

C.  FAULTY  COMPONENT  SELECTION 

This  program  was  written  to  determine  possible  and  probable  component 
failures  from  the  two  input  signatures.  The  signature  from  V,/I,  and 
Vp/V-,  are  both  processed.  A  component  is  flagged  as  possibly  failed  if 
one  of  the  library  signatures  for  that  component  is  selected  as  a  match. 
If  two  or  more  library  signatures  for  a  component  are  selected  as  a  match 
that  component  is  flagged  as  probably  failed. 

The  flow  graph  for  the  program  is  shown  in  Figures  22-24.  The  main 
program  allows  for  several  signatures  to  be  read  in  and  sent  to  subroutine 
DET,  one  pair  at  a  time. 

Subroutine  DET  first  checks  to  see  if  the  signature  indicates  the 
circuit  is  functioning  properly  (i.e.,  signature  is  all  5s).  If  so,  a 
message  is  printed  and  control  is  returned  to  the  main  program.  The  first 
entry  into  DET  with  a  non-nominal  signature  results  in  the  library  signa- 
tures being  read  in  and  stored. 

When  it  has  been  determined  that  non-nominal  signatures  are  being 
processed,  both  input  signatures  are  compared  with  the  library  signatures 
to  determine  the  closest  match  distance.  The  signatures  are  compared  a 
second  time  and  those  library  signatures  that  are  the  minimum  distance 
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Figure  22.  Flowgraph  of  Main  Program  for  DET 
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Figure  23.  Flowgraph  of  Subroutine  DET 
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Figure  24.  Flowgraph  of  Subroutine  DET  (Continued) 
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from  the  input  are  used  to  set  flags  for  determining  possible  and 
probable  faulty  components  as  described  earlier. 

The  program  is  listed  showing  the  nearest  neighbor  matching 
procedure. 

D.  SUBROUTINE  LSQ 

Subroutine  LSQ  was  written  to  determine  which  library  signatures 
most  closely  matched  the  input  signature.  LSQ  is  similar  in  instruction 
flow  to  subroutine  DET  and  therefore  a  flow  chart  is  not  presented. 

In  order  to  determine  a  match,  LSQ  uses  a  method  similar  to  the 
nearest  neighbor  concept.  However,  if  one  considers  the  signature  to 
be  a  vector  with  each  test  frequency  a  dimension,  LSQ  does  not  weigh 
every  dimension  equally.  If  one  dimension  can  be  strongly  affected  by 
random  variations  in  component  values  that  dimension  is  not  considered 
equally  with  the  other  dimensions.  Thus,  the  distance  between  the  input 
vector  and  the  library  vectors  becomes: 

,2 


,f,  ai(xi  -  V 


where 


t-bfy  +  b2x2  +  b3x3  +  ••• 


=  input  vector 

=  C1Y1  +  C2Y2  " 
=  library  vector 


Y=C1Y1   +C2Y2+C3Y3+   ••• 


a.  =  weighting  factor  for  the  i   dimension  (or 
frequency) 


a  •  <_  1  . 


5S 


LSQ  operates  on  only  one  signature  at  a  time  and  selects  the  proper 
library  signature  to  match  the  input.  The  output  is  a  list  of  the  pos- 
sible signatures^  and  the  component  variation  defined  by  each  one. 

One  may  choose  to  set  each  a.  =  1.  This  would  be  the  case  if  the 
quantization  process  was  sufficient  to  take  care  of  the  random  variation 
problem. 
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VIII.  SUGGESTED  FURTHER  STUDIES 

A.  ALTERNATE  APPROACH  TO  SIGNATURE  TABLE 

Table  VIII-1  lists  the  calculated  response  at  each  test  frequency 
for  the  primary  test  circuit.  Table  VI I 1-2  lists  the  range  of  the 
response  at  each  frequency  and  the  maximum  deviation  from  nominal  response, 
Zero  output  is  not  included. 

Test  Frequency  (Hz)  Response 

.00048 
.027 
.055 
.299 
2.99 
49.97 

TABLE  VIII-1.  Nominal  Response  at  Test  Frequencies 
Test  Frequency  (Hz)     Response  Range        Maximum  Deviation 


Vvi 

V1! 

-2.28 

Db 

12.86  Ohms 

-6.77 

Db 

3.22  Ohms 

-10.85 

Db 

14.32  Ohms 

r56.8 

Db 

1. 01   Ohms 

-120.3 

Db 

4.34  Ohms 

-215.5 

Db 

53.46  Ohms 

V2/V.,(db) 

V^I^n) 

VyV^db) 

yi^ 

.00048 

.0028 

-29.91 

3.0         66.77 

-27.63 

+53.91 

.027 

15.58 

-46.16 

.246  300.2 

-39.39 

+296.98 

.055 

-4.36 

-55.7 

3.0       302. 

'    -44.85 

+287.68 

.299 

-20.94 

-97.14 

.23     300. 

-40.34 

+296.99 

2.99 

-52.28 

-160.41 

3.         319.6 

+68.02 

+315.26 

49.97 

-98.48 

-255.5 

3.       5338. 

+117.02 

+5284.54 

TABLE  VIII-2.  Response  Range  and  Maximum  Deviation 

One  can  see  that  the  responses  do  not  cover  the  same  range  of  values  nor 
are  the  maximum  deviation  from  the  nominal  response  the  same.  From  the  nature 
of  the  network  functions  it  can  be  assumed  that  the  responses  at  different 
test  frequencies  are  not  independent  for  a  given  circuit  condition.  Due  to 

these  situations  it  would  seem  logical  to  weigh  equal  changes  at  different 

frequencies  differently. 
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The  setting  of  quantization  levels  and  subsequent  signature 
generation  is  one  way  to  accomplish  this  weighing.  The  setting  of 
quantization  levels  based  on  the  selected  situations  tends  to  weigh 
responses  at  some  test  frequencies  less  than  others.  However,,  this 
is  purely  an  "educated  guess"  method  and  has  no  analytical  foundation. 

An  alternate  scheme  was  studied  briefly  and  is  suggested  as  the 
subject  of  further  study. 

If  one  considers  each  test  frequency  as  one  of  the  dimensions  in 
an  n-dimensional  space  then  each  response  would  be  a  point  in  the  n- 
dimensional  space.  Instead  of  limiting  calculations  to  a  few  selected 
circuit  conditions  one  could  generate  many  points  representing  many 
circuit  conditions.  One  could  then  have  a  large  enough  sample  of  data 
to  do  some  statistical  operations. 

In  order  to  remove  dependence  the  normal  procedure  is  to  multiply 
each  vector,  representing  a  point,  by  the  inverse  of  the  co-variance 
matrix.  The  result  is  the  collection  of  points  in  a  "transformed  space" 
where  dependency  is  accounted  for. 

The  problem  is  the  construction  of  the  co-variance  matrix.  However, 
since  the  data  base  is  now  sufficiently  large, an  estimate  of  the  matrix 
can  be  calculated  from  the  data. 

With  a  circuit  under  test  the  response  would  constitute  a  vector 
which  would  be  multiplied  by  the  inverse  of  the  covariance  matrix.  At 
this  point  there  are  several  options  for  fault  identification.  One  can 
make  a  comparison  with  the  stored  data  and  make  the  fault  identification 
using  the  nearest  neighbor  technique.  An  alternate  approach  would  be  to 
use  an  average  of  the  nearest  m  neighbors  when  m  is  some  positive  integer. 
There  are  other  possible  procedures  for  matching. 
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If  one  is  only  interested  in  which  component  is  faulty  it  might 
also  be  of  interest  to  investigate  any  possible  grouping  of  points  which 
would  be  useful  as  in  pattern  recognition.  The  unknown  could  then  be 
assigned  to  a  group  according  to  its  position.  Reference  5  details  the 
necessary  pattern  recognition  information. 

The  above  described  procedure  would  require  a  computer  with  sufficient 
storage  to  process  and  store  the  information  pertaining  to  the  data  base 
and  covariance  matrix.  However  such  a  computer  is  readily  available  for 
a  nominal  cost. 

B.  EXTENSION  OF  PROCEDURE 

It  soon  becomes  obvious  that  for  even  moderate  sized  circuits  the 
procedure  would  be  wery   difficult  to  use.  For  example,  a  25-element 
circuit  with  6  possible  element  value  perturbations  would  yield  150 
signatures.  For  ease  in  testing  the  number  of  test  frequencies  should 
be  as  small  as  possible.  However,  with  150  signatures  considerable 
signature  duplication  is  probable  unless  the  number  of  test  frequencies 
is  large. 

Reference  3  shows  how  a  ladder  network  may  be  divided  into  subnet- 
works such  that  the  voltage  gain  of  a  given  subnetwork  is  independent  of 
the  voltage  gain  of  all  preceding  networks.  Therefore,  testing  by  add- 
ing subnetworks  one  at  a  time  can  yield  the  fault  component. 

This  procedure  uses  only  the  voltage  gain  but  does  show  one  approach 
to  handing  large  circuits.  In  any  case,  a  large  network  would  have  to 
be  divided  into  subnetworks  by  some  scheme  in  order  to  reduce  the  problem 
to  a  manageable  size. 
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IX.  CONCLUSIONS 

The  majority  of  the  conclusions  made  as  a  result  of  this  study  were 
based  on  a  comparison  of  the  results  is  section  VI  and  results  published 
by  other  investigators. 

Most  investigators  reported  approximately  75%  isolation  using  one 
network  function  with  no  allowance  for  van* ation-wi thin-tolerance  of  the 
nonfaulty  elements.  One  investigator  [2]  reported  78%-98%  success.  How- 
ever, the  results  were  not  sufficiently  documented  to  judge  the  merit  of 
the  claim. 

The  results  of  this  study  show  that  an  additional  error  rate  of 
17%-42%  can  be  introduced  by  the  within  tolerance  variations  if  the 
conventional  procedures  are  followed.  However,  if  sensitivity  functions 
are  used  to  aid  in  selecting  test  frequencies  the  error  rate  due  to  other 
element  variations  drops  to  11%-31%  for  one  network  function.  Further- 
more, using  two  network  functions  fault  isolation  is  successful  for  up 
to  94%  of  random  single  fault  conditions. 

The  high  degree  of  successful  fault  isolation  clearly  points  out  the 
advisability  of  using  at  least  two  network  functions  for  testing.  The 
two  functions  used  for  this  study  were  V?/V\  and  V-./I, .  However,  the 
designed  purpose  of  the  circuit  may  dictate  a  different  choice  of 
functions  to  be  used. 

The  results  clearly  show  the  improved  accuracy  of  isolation  using 
sensitivity  functions  as  a  guide  in  test  frequency  selection.  In  addition 
to  improved  isolation  this  method  also  eliminates  the  possibility  of 
unknowingly  picking  test  frequencies  which  yield  the  same  information. 


63 


COMPUTER  OUTPUT 
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Figure  25.  Signatures  from  Basic  Procedure 
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Figure  26.  Signatures  using  Sensitivity  Procedure 
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COMPUTER  PROGRAM  LISTINGS 

REAL**  UMcp,A(320) 

RFAL  MAG,NMAG,WMAG(2»320) 

DIMENSION  p«cr,n  (  5  ,  5  )  ,  c(  1  0  )  ,D(  ?  )  ,0  (  5  )  ,.R  (  c)  ,  F  (  5)  ,  HN(  5  ) , 
1Y(5IiEI(5)  ,011  10),  MAG  (4, 10)  »  NMAG  (  6  ,  10  ,  ^0  )  ,r»x(5),YX(5), 
2Z(  5)  , X"AX(  20),XMIN(20),INMAG(6,10,20)  ,,OMAG<  2»10) 

COMVPM  NMar; 

DATA    OMFGA/?c*0.  /,  0/^*0.  /,R/5*0./,  F  17*5*0../  ,F/ 5*0./ 

DATA    01/10*0./ 

READ    (5,110)     P4CQ 

F( 1)=3. 

F( 2)  =  . 17 

E( 3)=1  .0000^1 

E(4)=10. 00001 

F(5)=9. 

F(6)=4. 

DO    3*    1=1,6 

DO    3P    K=l,10 

GO    TO    ( 1,2, 3, 4, 5, 6, 7, 8, 9,10), K 

1  TFMP=C( I ) 
GO    TO    11 

2  F(  I )  =  1  ,1*TFMP 
GO    to    11 

3  F(  I  )=0.9*TEVP 
GO    TO    11 

4  F( I )=0. 
GO    TO     11 

5  E(  I  )=0  .C*TC,-'D 
GH    Tp     11 

6  F(  I  )=0.3*TEMP 
GO    TO    11 

7  F(  I)  =  l  .2*""«^P 
GO    TO    11 

8  F(  I )=1  ,5*TFMp 
GO    T0    n 

9  F( I  )  =  100.*TFMO 
GO    TH    1  1 

10  F(  I  )=TFvr> 
GO    TO    3  9 

11  UN  < 1 )  =  F  C  4 ) 

D< 1 )  =  F( 1 )+F( 4) 

D(2)  =  F(2)+F(c5)+r<l)*c(M*(F(3)4-E(5)) 

D(3)=c(5)*(-(l)*-(3)+F(M*-(6))+c(4)*E(7)*(E(3)+F(5)) 

D(5)=E(2)*E(  3)*F(4)*t=(5)*E(6) 

Y( 1)=1  . 

Y( 2)=c( 4)*(c(3) *c( M ) 

Y(^)=c( 3)*c( c) 

Y(4)=E(3)*E( ^)*F(5)*F(6) 

IF     (K.GT. LOR. T.OT.l)     OH    TO    36 
CO  WRITF     (8,100)    UM(1) 

CC  WRITr     (B,10l ) 

CC  WRITF     (8,102)     (0(N),N=1,5) 

CC  WRITF     (P,10?)     (D(M),N=1,5) 

CC  WRItc     (Qtirv. ) 

CC  WRIT-     (8,103)     (Y(N),N=1,4) 

C    DFLFTF    STARTING    HERF    FOR    SENSITIVITY    FUNCTION    USE 

00    1111    LM=1 ,5 

DX(LM)=D(LM) 
1111    YX(L*M=Y(LM) 

CALL     PROD    (nx,c,otEI  ,PCL,  IRf  I«=R)    • 

IP    (  IFR.rO.O)    GO    TO    12 
CC  WRITF     (8,105)     I-c 

CC    12    WRITE     (*,106)     (Q(N),EI <N) ,N=1,4) 

12  CONTINUE 

CALL     P»OP    (YX,4,d ,,FtpnM, ipt IER) 

IF    (  IER  .cQ.O)    on    Tn    13 
CC  WRTTr.     (?,10«-)     Tco 

CC     13    WRITF     (8,106)     (o (N) ,F(N) ,N=1 ,3) 

13  CONTINUF 

DO    15    M=l,A 

IF    ( &BS( BT(M) ) .LF. 1=-10)    GO    TO    14 

OMEGA     (1,")=ABS(FI(M) ) 
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14 

15 


16 

17 


75 

ib 


19 


20 
21 


22 
23 


24 


25 

26 


27 
2* 


GO  TO 
OMEGA 
CONTIN 
00  I7 
IF  (  AB 
OMEGA 
GO  th 
OMEGA 
CONTJN 
DO  IB 
00  1^ 
00  IP 
IF(OME 
GO  TO 
OMEG& 
CONTIN 
DO  21 
DO  2  0 
IF  (OM 
GO  TO 
OMEGA 
OMEGA 
OMEGA 

CONTTNJ 

CONTIN 

DO  2  2 

IF  (OM 

LP=L+1 

OMEGA 

OMEGA 

OMEGA 

IF  (  nv 

LP2=L+ 

IF  (LP 

OMEGA 

OMFGA 

OMEGA 

IF  (LP 

LP3=L+ 

OMFGA 

OMEGA 

OMEGA 

GO  TO 

CONTIN 

DO  2  6 

00  25 

IF  (OM 

GO  TO 

OMEGA 

OMEGA 

OMEGA 

CONTIM 

CONTIN 

DO  2  7 

IF  (  DM 

LP=L+1 

OMFGA 

OMEGA 

OMEGA 

OMEGA 

IF  (C^ 

LP2=L+ 

IF  (LP 

OMEGA 

OMFGA 

OMEGA 

gr  th 

CON^IM 

M2  =  Q 
DO  2  9 
DO  29 

M2=M2+ 


15 

( 1  ,M)=A8S(0(M) ) 

OF 

S( 
(2 


1.3 

F(  v) >.LF.lc-io 

,m)=ABS< F( M) ) 


)  GO  TO  16 


(2,M)=ABS(R(M) ) 

UE 

M- 


NO 
GA 
IB 

(M 

OF 
J  = 
L  = 
EG 
20 
(! 
(1 
(1 
OF 
Uc 
L  = 
FG 


1,2 
It  4 
=  1,4 

( ^,N) .EG. OMEGA (M,NO) .AND.N.NF.NOJ  GO  TO  75 

,N)=0. 

1,4 

J, 4 

A(l,  J)  .GT.OMFGAn  ,L)  )  GO  TO  19 

,5)=0MEGA  (l»J) 
,J)  =  OM=GA  U,L) 
tL)=OMEGA  (1,5) 


1,4 

A( 1,L) .CC0. ) 


(3, 

(3, 

(^, 

EGA 

7 

2.G 

(  3, 

(3, 

(3, 

2.E 

3 

(3, 

(3, 

(3, 

23 

UE 

J  =  l 

L  =  J 

EGA 

">*> 

(2, 
(2, 

(2, 
UE 

Uc 
1=1 

EGA 

(^, 
(4, 
(4, 
(  *, 

EGA 

2 

2.G 

(4, 

(4, 

(4, 

?° 

UE 


r>  )=(OvcGA  (  ltL 

1  )=OMEGA  d,L 

2)=0MFGA  (1,L 

(3,1 ) .LT.O. ) 

T.4)  GO  Tn  yx 
c  j-(HMPr,A(  l  ,C 
3  IOMEGA  (  1  ,  L 
^)=OMEGA  C1,L 
0.4)  GO  TO  23 

S)=(OVPGA  (1  , 

4)  =  0MEGA  (1,1 

5)  =  nvcr7A(  i  ,  LP 


GO  to  22 
P)_nMPGMl,L)  J  /2. 

1-nMcr,*   (">?^) 
)+rrMcr,A  (3,5) 
hmcGA  (3,l)=H«E0A(l,L)/2, 


p?)-^mCGA(1  ,LPH  S2. 
P)+ny  =  GA  (?,<=) 
t>2)  +  0MEGA<3,51 


LP3)-0^EGA  (l,LP2)jr2 
31+OMEGA  (3,5) 


,3 

,3 


(2,J).GT.0MEGA  (2,D)  GO  to  24 


5)  =  0MEGA  (2, J 
J)=OMEGA  (2,L 
L)=OMEGA  (2,5 


(2,L).E0.0.)  GO  TO  2n 


T.^)  GO  T^  ?* 
5)=(nMFGA  (2, 
3J  =  0M=GA  C,L 
4)=0MEGA  (2»L 


LP)_rvFnA  (2,L))/2. 

)-omcga  (4,5) 

)4-0MEGA  (4,5) 

oj+qv^A  (4,c?| 

OMEGA  (4,l)=OVEG4f2,L*/2 


LP2)-0MEGA  f?,t°n^2. 
D)+0"E0A  (4,') 
P2)+nvcGA  (4,5) 


M=lf4 

Ml=3,4 

1 


29 


291 


31 
32 


33 

34 

35 

C 

36 


37 


cc 

cc 
cc 


372 


4<50 
500 


371 
38 
39 

40 
41 

411 


CC  42 
42 

CC 

43 


CC  44 
44 


01  CM? ) 

01(^2+ 
DO  2  9! 
00  2Q! 
IF  <ni 
0X=01( 
IP  <nx 
CONTIN 
00  3  2 

00  31 
Ic  (ni 
OK  10) 

01  <M3) 
ni  (Mi) 

CONTTN 
OONTTN 
00  ?4 
IF  (01 
M6*10- 
IF     (M5 

00  3  3 

01  (^7) 
GO    to 
CO NT  IN 
0 1  (  V  6  ♦ 

■    STic 
CALL     w 
IF     (K. 

00    37 
MAG    (3 
MAG    (4 

WR  I  t  E 
WRITr 

WRT  TC 

WP  IT' 
PO  ^00 
DO  500 
IF(MAG 
OMAG(  I 
GO  TO 
DMAG(  I 
CONTIN 

00  38 
IF     (MA 
NMAG( I 

1  F  (  M  A  G 
NMAGf  i 

CONTIN 
00    3  9 
CALL    N 
WR  I  TF 

00  4  0 
CALL     0 
DO    4  1 
CALL     O 
DO    411 

nn  4ii 
DO    411 

1  N  M  A  G  ( 
DO    4  3 
nn   4? 

WR  I T  F 
CONKM 
WB  ITT 

CON  "KM 
DO    4  5 
DO    44 


sHMCRJ ( Ml  f  M) 

1  )=0MEGA(3,5) 
IL=ltlO 
IJ  =  T  L ,10 

(  I  J  J  .FQ„0. )     GO    TH    291 

m/nidj) 

.Mc. 1. .AM0.(nx.GT..Q.ANn.OX.LT.l.l) )     OK  IL)=0. 

M3=l,« 

(  "3)  .LP. 01  (V<V)  )    GO    to    31 
=  r i  ( m?  ) 
=  0 1  (  M  4  ) 

=01(10) 

uc 
uc 

M5=l ,9 

(M5) .r0.0. )    GO    TO    3* 

.50.1)    GO    TH    35 

M"»=l  f  MA 

=  0]   ( MS+M7-1 ) 

35" 

UF 

1)=0. 

DFLpTinN  nr-oc 

AGl    <MA,ni  ,UN, n,  y,  *AG  ) 

NS.l.OR. I. MP. I)     GO    TO    372 

L"=1,'0 

,LM)=MAG(  1,L«I 

,LM)=MAG(2fLM) 

(  4  , 1 1  5  )     M  *  r,  (  3  ,  L  M  )  ♦  MAG  (  4  ,  L  m  ) 

{  8,1  0"7)     (ni(  MR)  f  Mq  =  l  ,M6) 

<4,im     (  OK  mb)  ,M8=1  ,  V6) 

( 8»108) 

(8,10?)     (Mind  ,N),  MAG(2fN)  ,N  =  1  ,M6) 

IK=1,2 

IL=1,V6 
(  I*,  I!)  .Lc.0.  )    GO    TO    4QQ 
K,IL  )  =  20.*AL0010(MAG(  IK,  ID) 
500 

K,IL)=0. 
UF 

L"  =  l  ,10 

G( 3,LM) .LF.O.  )    GT    T0    371 
,  K  ,  L  M)  =m  *g  ( i ,  lm )  /MAG  (  3 , L M ) 
(4,LM)  .LF.O.  )    r^n    in    33 
,K,LM+10)=MAG(2tLM)/MAG(4,LM) 
U- 

K  =  l,l<3 
0*>M6,?,3,K,XM£X(f)  ,XMIM(K)) 
(4,H5)     (XMAX(K7.)fXMIN(KZ),K7=lf19) 
K  =  T      0 

UANT4     (K,XMIN(K) ,XMAX(K) ,PACO) 
K  =  1 1 , 1  9 
UANTA     (K ,XMIN(  K)  , X  VAX ( K )  , F AC  0) 

T=l,6 

J  =  l  ,10 

K  =  1  ,  ■>  0 
T  ,  J,K)=NMAG(  I  , J,K) 
1=1,6 
.)  =  <,, c 

(a,ll2)     (  INMAGU  ,.J,K)  ,K=1  ,M6) 
U^ 


(  p  ,  1 1 3  ) 
Uc 

1=1,6 
J=4,9 

+  10 

WRITE     (8,112)     (  IN^AG( I , J,K)  ,K  =  11  ,M6P) 
CONTINUE 
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cc 


45 


100 
101 

102 


WRIT 
CONT 
WR  IT 
STOP 

FOR" 


P  (  p 

INUF 

n  (  4 

1 


,113) 

,112)  (  (  (  INMAGM  »J,K)  ,K=1, M61,  J=^ta)  rT=l,6) 


4T  { 
FORMAT  ( 


1 5  X , 

6X,' 


1PC13.5) 
V2/V1  = 


103 

104 

105 
106 

107 
10« 
10° 
110 
11! 
112 
113 
115 
116 


P  0  R  M 
IX,' 

tlPF 

FORM 

**2 
FQpv 


AT      ( 

13  .c 
AT    ( 

+     '1 
AT    ( 


1  ) 

15X,lPE13.5flXf  «♦•  ,lPE13.5flX,  •$**!    +-»'IPE 
+    • 

,  lx, 
15X, 

pr  1  3 

6X,  • 


•S**3    ♦    «  ,IPE13.5  ,1  X,  «S***  »,./) 
1Pr!3.5,»    +    ,lPE13,5fl     S**T     +    rIPc:T3. 


Zl  .1    = 


FORM 
FORM 
RT     R 

FOOM 

CQQM 

c0RM 

FORM 
FORM 
FIRM 

FORM 

POP" 
pnpM 

STOO 

END 


AT 

AT 

FAL 

AT 

AT 

AT 

AT 
AT 
AT 

A  T 

AT 

4T 


( 
( 

( 

( 

( 

(  F 

( 

( 

{ 

( 

( 


110) 
//// 
PART 
//// 

//:  x 

1PE1 
F5.0 
IX, 1 
IX,7 
1H  ) 
?F1A 
5F15 


) 


10X,,Rc^L     PA^T     IMAG    PART    RFAL     PART    TM 
IMAG     P&PT     Rh^l     PAPT     TMJG     PtOTl/QX,flP 
IX, 'THE    FOF^UFMCIFS    USED    ARF: ♦ t / 1 °} oc 
,M/?/Vl     MAGNITUDE  Zltl        MAGNITUD 

c..c,  1PEP1.5) 
tlc) 

P020.5,1P2E20.51 
T4) 

.?) 
.7) 


1 3., 5, 
5,."  S 


AO  DA 
10., 3) 

E«  /) 


SUBROUTINE  MAGI  (  N  ,  OMEGA  ,  XNUM,  DEN, YN'JM,»*AG  ) 
RFAL  mag 
D I MENS  10 
0^  1  1=1 

w=omfga< 

YP  =  YNUM{ 
Yl=YNUM< 

OP=DFN(  l 

DI=DEN<  2 

YM  =  SOR  T{ 
. — GM  =  SOPT{ 

MAG( 1,1) 
1    MAG(?,T) 

RETURN 

END 


\|    OMEGA (10)  ,XNUV(5),DEN(5)  ,YNUMfc)  tMAG<4, 10) 

,N 

I  ) 

I  )-YM(JM{  1)  ftM*w 

->  )  *  w-  Y  N I  i M  (  &  )  *  W *  *3 

)-OcN(  3)  *W*W4-0FN(  5)*W**4 

) #w_0FN( 4) *W**3 

YP  *YR+YI  "Y  T  ) 

r)p  *np+n t *ot  ) 

=  XNU*M I ) /OM 

=nw/ ym 


SUBROUTINE    QUANTA     ( K , OL, OH , F AC ) 
C        K=TFRM     IM     ARRAY 
C       0L=10WFR    LTmtt    poo    trl. 
C       QH=UPPFR    LIMTT    FOR     tol. 
C       FAC=FACTnQ    cod    QUANTIZATION 

REAl     MWArT(6,iot20»  ,Z(  5) 

COMMON    NMAG 

no    5    1=1,6 

DO    5     J =1,10 

QU-NMAG< I  ,J,K) 

7( 1)=PAC 

7( 2)=FAC*3. 

7( 3)=FAC*6. 

Z(4)=FAC*12. 

IF     (O(l.LT.OL)     GO    TO    3 

IF     (OU.GT.OH)     GO    TO    1 

NMAG< I ,J,K)=0. 

GO    m    5 
1    DO    2    L  =  1 ,  A 

OF AC =7 (L )*( OH- I. )+0M 

IF     (OU.GT.QFAC)     GO    TO    2 
J,* )=L 

J,K)=5. 


NMAG 

(  1* 

GO 

TO 

«5 

NM, 

AG 

(  T, 

GO 

TO 

5 
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3  DO  4  M=l,« 

OFAC=OL-7(w)*( l.-Ol ) 
IF  (QU.LT.OcAC)  HO  TO 
N"AG  ( I , J, K) =M+4 
GO  TO  5 

A  NMAG( I »J,K)=10. 

5  CONTINUE 
RETURN 
END 


SUBROUTINE    MOM 
PEAl     NMAG(6»10f 
COMmhn    mm AG 
X«AX  =  1  . 
XMIN=1  . 
DO     1    L  =1  ,  I 
DO     1     M=J1,J? 
IF    (NMAG(LvMyKl 
IF    (NMAG(L,M,K) 
CONTIMUc 
RETURN 
END 


(  I  ,  Jl 

20) 


J2tK,XMAX,XMIN) 


GT.XMAX)     XM^X=N^AG(L*M-rK) 
LT.XMIN)     XMlN=NMAG(L».Mf.K) 
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DIMENSION    ISIS(l2/9/6)/ JSI3C6) 
DO    1     1=1/6 
DO    1    K=i/9 

1  READ  (5/ 100)     (  ISIS(  I/KiU>Lsl*6)#(  rsT3(  I>6/;K-/.>D  /  L*  \>  6  ) 
DO    <*    !«1/12 

DO    <*    U=l/9 
09    3    K*l/6 

3  JSI3(K)  =  ISI3(  I/J/O 

WRITE (6/ 101 )  I/J/(USlG(Kl)/Kl=l/6) 
CALL  LSQ(I#Ji JSI3) 
k    C0NTINUE 

100  F9RMAT( 1212) 

101  F0R*AT(//5X/  'THF  INPUT  JS:  ELEMENT**  £3*''   VAR.- •  N9i  » /  I  2 
END 

SUBROUTINE  LS3 ( I  1 / Jl / US  I  3 ) 

REAL  ID(36)iIS 

DIMENSION  ISI3(36#6)i  JSIGC6>#<3IG(36*6)/.-FAC(6l-) 

IF  {  11-1)  1,6*1 

6  IF  (Ul-1)  1/9/1 
9  FAC(l)»lt 

FAC(2)«.3333 

FAC(2)*.3333;FAC(3)  =  .565;FACU)=*^62;,FAC(5)-«<69> 

FAC(6)*«7^5 

DO  2  J=l/36 

2  READ  (5/ IOC)  (ISIG(J/I)#I«l#fr)!*tKStG(;j>r)'//IM^'6-): 
1  IF(H-7)  3/7/7 

7  DO  5  U*l/36 
ID(J)-0 

DO  k    1=1/6 

1D1»KSIG<  J/  I  )-USH(  I  ) 

4  ID(J)«IO(J)+FAC(I )*ID1**2 

5  CONTINUE 
GO    TO    17 

3  DO    15    U=l/36 
ID(J)*0 

DO    10    1=1/6 

1D1MSIG(U/I  )-JSl3(  I  ) 
10    ID(U)=ID(U)+FAC( I )*I01**2 
15    CONTINUE 
17    IS=ID(1) 

DO    20    1=2/36 

IFdD(J)-IS)    21/21/25 

21  IF {  1 1-7  )    22/8/8 

8  WRITE(6/101)    Ji  (K5IGU/  U)#TJ«l/6) 
GO    TO    25 

22  WRITE(6/101)  J/ (  ISIG(J/IJ)#IJ  =  U6) 
25  CONTINUF 

100  F9RMAT(6I2/ 16/512) 

101  F0RMAT(5X/'A  PSSSI3LE  SIGNATURE  lStv»Z3ts>5Xi&l'Z.Y 

102  FOR^AT( 1H  ) 
WRITE(6/102) 
RETURN 

END 
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MAIN  PROGRAM  FOR  USE  WITH  SUBROUTINE  DET    ( FAC  NOT  INCLU 
DIMENSION  ISIG(12,9,6),JSIG(6),LSIG(6) 
DO  1  1=1,6 
DO  1  K=l,9 
1  READ(2,100)  ( ISIG(I,K,L)  ,L  =  1,6), <  I  S  IGU +6,  K,  L  )  ,L  =  1 ,  6) 
DO  4  1=1,6 
DO  4  J=l,9 
DO  3  K=l,6 
JSIG(K)=ISIG(I,J,K) 

3  LSIG(K)=ISIG( 1+6, J, K) 

WRITE(3, 101)  ( JSIG(K) ,K=1,6) , (LSIG(K) ,K=1,6) 

4  CALL  DET( I , J , JSI G, L SIG ) 
STOP 

100  FORMAT(12I2) 

101  F0RMAT(/,5X, 'THE  INPUTS  ARE : ' 1 4 , 5  I  2 , •  -  -  ',612) 
END 

SUBROUTINE  CET(  I,  J,  JSICLSIG) 

DIMENSION  ISIG(36,6),KSIG(36,6),JSIG(6),LSIG(6),ID(36) 

DO  3  K=l,6 

IF( JSIG(K)-5)  13,12,13 

12  IF(LSIG(K)-5)  13,3,13 

3  CONTINUE 
WRITE(3,101) 
RETURN 

13  IF(I-l)  2,23,2 
23  IF(J-l)  2,33,2 
33  DO  1  K=l,36 

1  READ(2,100)  ( ISIG(K,L),L=1,6) , (KSIG(K,L) ,L=1,6) 

2  DO  4  K=l,3  6 
ID(K)=0 
LD(K)=0 

DO  4  L=l,6 

ID1=ISIG(K,L)-JSIG(L) 
ID2=KSIG(K,L)-LSIG(L) 
ID(K)=ID(K)+ID1**2 

4  LD(K)=LD(K)+ID2**2 
DO  5  L=l,6 

5  N(L)=0 
IS=ID(1) 
IT=LD( 1) 

DO  7  11=2,36 

IF(LD( II )-IT)  16,17,17 

16  IT=LD(I1) 

17  IF( ID( II )-IS)  6,7,7 

6  IS=ID(I1) 

7  CONTINUE 

DO  10  K=l,36 
IX=(K-l)/6+l 
IFdD(K)-IS)  35,35,36 

35  N( IX)=N( IXJ+1 

36  IFCLD(K)-IT)  37,37,10 

37  N( IX)=N( IX ) +1 
10  CONTINUE 

DO  15  K=l,6 
IF(N(K)-2)  40,45,45 
40  IF(N(K)-1)  15,46,15 

45  WRITE(3,102)  K,N(K) 
GO  TO  15 

46  WRITE(3, 103)  K 
15  CONTINUE 

RETURN 

100  F0RMAT(6I2, 16,512) 

102  FORMAT(«     THE  FAILURE  IS  PROBABLY  COMPONENT » , I  3 , 1 1 0 ) 

103  FORMATC     THE  FAILURE  IS  POSSIBLY  COMPONENT ',  1 3  ) 

101  FORMATS     THE  SIGNATURE  IS  NOMINAL') 
END 
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